Quantum Spin Hall, triplet Superconductor, and topological liquids 

on the honeycomb lattice 



Cenke Xu 

Department of Physics, University of California, Santa Barbara, CA 93106 
(Dated: October 5, 2010) 

We classify the order parameters on the honeycomb lattice using the SO (4) symmetry of the 
Hubbard model. We will focus on the topologically nontrivial quantum spin Hall order and spin 
triplet superconductor, which together belong to the (3,3) representation of the SO(4). Depending 
on the microscopic parameters, this (3, 3) order parameter has two types of ground states with 
different symmetries: type A, the ground state manifold is [S 2 <S) S 2 \/Z2\ and type B, with ground 
state manifold SO(3) ® Zi. We demonstrate that phase A is adjacent to a Z2 ® Zi topological 
phase with mutual semion statistics between spin and charge excitations, while phase B is adjacent 
to a nonabelian phase described by SU(2) Chern-Simons theory. Connections of our study to the 
recent quantum Monte Carlo simulation on the Hubbard model on the honeycomb lattice will also 
be discussed. 



I. INTRODUCTION AND SYMMETRY 

We consider a class of (extended) Hubbard models on 
a bipartite lattice at half filling, with the following form: 

H = -t4,a c j,<r+H.c. + Un itt n u +H' 

<i,j>,a 

h' = %4^> + H - c - 

iGsA, jEsB 

~t" ^ JijSi • Sj -\- VijTi • Tj -\- 

i,j 

% = ((-l) l Re[A] l7 (-l) l Im[A],. m - l) , (1) 

Si = ^c\(7Cj is the spin operator, A, = c\ia v Ci is the 
on-site spin singlet Cooper pair. This extended Hubbard 
model has a manifest SU(2) spin symmetry. However, 
if electrons only hop between two different sublattices 
(sublattice A and B, denoted as sA and sB), after a 
sublattice depended particle-hole transformation for spin 
down electrons: 

cu -> (-l)'cL' (2) 

the Hamiltonian is almost unchanged except that U 
changes sign, and Jy switches with V%j. This implies 
that in addition to the apparent SU(2) sp i n symmetry, 
this model also has an SU(2) c h ar ge symmetry that mixes 
(ci.-|-, (— l)* c !^)- Therefore the full symmetry of this ex- 
tended Hubbard model is 1 ' 2 

SO(4) ~ [SU(2) spin ® SU(2) chargc ]/Z 2 , (3) 

for arbitrary parameters in Eq. [TJ For instance, this 
SO (4) symmetry holds for the simplest Hubbard model 
with only on-site Hubbard interaction and nearest neigh- 
bor electron hopping. 

Since the SO (4) symmetry is the full symmetry of the 
extended Hubbard model Eq. [ljon any bipartite lattice, 
all the order parameters should be classified in terms of 



the representations of the SO(4) Lie algebra. In this 
work we will take the honeycomb lattice as an example. 
Using the notation introduced in Ref£, we expand the 
electron at two Dirac valleys by d\^ — e z< ^ ia ' r c (where 
Q1.2 = ^(^TJjO) are the wavevectors of the valleys), 
and introduce Pauli matrices t° and fi a which act on 
the sublattice and valley spaces respectively. Then, after 
introducing real Majorana fermions £ a as the real and 
imaginary parts of e l ^ T e 1 ^^ (di, ir^d-if we obtain the 
continuum Lagrangian for the semimetal phase 

8 

Co=J2^M«- ( 4 ) 

a=l 

Here /1 is a 2+1 dimensional spacetime index, and the 
Dirac 7 matrices are (70,71,72) = t v ,t z ,t x , ( — C*7°- 
Using this notation, in the low energy field theory, the 
SU(2) S pi n and SU(2) c har g c symmetries are generated by 
the following matrices^: 

S x = o- x p v , S v = cr v , S z = a z p v , 

T x = a v p z , T v = a v P x , T z = p y . (5) 

a a are spin Pauli matrices, while p a are Pauli matrices 
that mix the real and imaginary parts of electron. Notice 
that SU(2) <g) SU(2) is a double covering of SO (4), which 
leads to the Z2 in Eq. [3] 

Based on the symmetry Eq. [3] and Lie Algebra Eq. [5j 
the spin and charge are dual to each other for a large 
class of extended Hubbard model. This spin-charge du- 
ality will lead to many interesting results in our paper. 
The structure of this paper is as following: in section 
HT1 we show that the two types of topological orders, the 
quantum spin Hall (QSH) order and triplet supercon- 
ductor (T-SC) are unified as one representation of the 
SO (4) group, and the Ginzburg-Landau theory gives two 
types of ground states with different symmetry breakings. 
Section III and IV study the phase diagrams driven by 
proliferating the topological defects in the two types of 
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between next nearest neighbor sites: 




FIG. 1: Honeycomb lattice and the vectors e a . 



orders described in section II respectively. In both sec- 
tion III and IV, we will first give an argument of the 
phase diagram based on the quantum numbers of the 
topological defects, then a more solid description based 
on the Majorana liquid formalism developed in Ref£ will 
be presented, and the results from these two approaches 
match perfectly with each other. Section V discusses 
the situation with SU(2) c h arg e broken down to U(l) c h a rge 
symmetry. 



II. SO(4) CLASSIFICATION AND 
GINZBURG-LANDAU FORMALISM 

Using the SO (4) algebra in Eq. [SJ we will classify the 
order parameters which immediately open up a mass gap 
for the Dirac fermion in the semi-metal phase. Some 
simple Dirac mass gap order parameters can be classified 
with these symmetries straightforwardly. For instance, 
the quantum Hall order parameter is a (1,1) repre- 
sentation of SO(4) i.e. it is a singlet of both SU(2) sym- 
metries. The two sublattice Neel order N a = £S a [i y ( 
is a (3, 1) representation. The fermion bilinear M a = 
(T a p y ( which belongs to the (1,3) representation is a 
two sublattice charge-density wave (CDW) and s— wave 
superconductor: M z ~ (-l)*(n< - 1), M x + iM y ~ 
c\icr v Ci. M a can be viewed as the spin-charge dual ver- 
sion of N a . In the simplest Hubbard model, iV a and 
M a orders can be realized in the two limits U \t\ and 
—U^> \t\ respectively. 

Now we discuss the following order parameters which 
belong to the (3, 3) representation of SO(4): 



lab — (A a bC, 



A ab = T a S b = 



/ —a z p x p z a x p x \ 

a z p z p x —a x p z 
V <t x a v p y a z J 



(6) 



This 3x3 matrix Q a b has SO(3)i c ft and SO(3) r i g ht 
transformations, which correspond to SU(2) c h a rge and 
SU(2) S pi n symmetry respectively. Qzb corresponds to the 
QSH vector^, while Q 2 b + iQ\b is the spin triplet pairing 



hb ~ £ 

jesA, a=l,2,3 

- (sA -> sB) 



Q2b + iQib ~ £ ic t J ii7 v a b c j+ea 

j£sA, o=l,2,3 

+ (sA -> sB). 



(7) 



ei = \/3x, e2, e 3 = — ± |y are three vectors on the 
honeycomb lattice that connect next nearest neighbor 
sites. Therefore the two types of topological order pa- 
rameters, QSH and T-SC, are unified through the SO (4) 
symmetry. Under time reversal symmetry T, and 
Q2b are even, while Qu is odd. Under reflection symme- 
try P x : y — » —y, Qu, and Q26 are even, while Q^b is odd; 
under P y : x — > —x, all components of Q a b are odd. 

The low energy dynamics of Q a b can be described by 
the following Ginzburg-Landau field theory: 



L Q = tr[a /i Q'9 / ,Q]+r(tr[Q t Q])+.g(tr[Q'Q]) 2 
+ u(ti[Q t QQ t Q]) + --- 



(8) 



The first three terms have an enlarged SO (9) symmetry 
which corresponds to the rotation between the nine order 
parameters in matrix Q; while the last term ~ u breaks 
this SO(9) symmetry down to SO(4) symmetry. Another 
term Det[Q] is also invariant under SO (4) transforma- 
tion, but Det[Q] breaks the time- reversal and reflection 
symmetry of the honeycomb lattice, therefore Det[Q] is 
forbidden in the Lagrangian Eq. [51 However, if the sys- 
tem already breaks the time reversal and reflection sym- 
metry (for instance (££) 7^ 0), Det[Q] would be allowed. 

In Eq. [8j when r < 0, Q is ordered, and the SO(4) 
symmetry is broken down to its subgroups. Depending 
on the sign of u, there are two types of ground states: 

Type A, u < 0, one example state of this phase is 
{Q33) 0, and all the other components (Q a b) = 0. In 
this phase the SO (4) symmetry is broken down to its 
following subgroup: 



( U(l) spin <g> U(l) ch argo ® Z 2 ) /Z 2 



(9) 



The U(l) sp in and U(l) c hargo symmetry are generated by 
matrices S z and T z in Eq. [5] The Z 2 in the numera- 
tor corresponds to reversing the direction of S z and T z 
simultaneously, while keeping Q33 invariant. The Z 2 in 
the denominator is the same Z 2 as in Eq. [3j which cor- 
responds to changing the sign of electron operator. The 
ground state manifold (GSM) of this phase is 



GSM ~ [5 spin ® Charge] I %2 



(10) 



The ground state can be described by two independent 
unit vectors N s and N c which belong to the (3, 1) and 
(1,3) representation of SO(4) respectively, and Q a b — 
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N^N b . This phase has four independent Goldstone 
modes. The Z 2 in Eq. [TU]is due to the fact that N s and 
N c can reverse direction simultaneously, and the ground 
state remains invariant. 

Type B, u > 0, one example state of this phase is 
(On) = (Q22) = (Q33) ¥= 0. The SO(4) group element 
G so(4 ) can be written as G su ( 2 ) :Sp in <8> G su ( 2 ), charge) and 
the type B phase breaks the SO(4) down to its subgroup 
with G su(2 ), S p in = ±G su (2) :Ch argc- This implies that the 
residual symmetry group elements can be parametrized 
as ±R a (9 a ,n s ) with 9 s £ (0,2n), which is equivalent to 
the diagonal subgroup SU(2) + generated by operators 
G a = S a + T a . Here R s (8 s ,n s ) represents spin rotation 
by angle 9 s about axis n s . The GSM of phase B is 

GSM = ( SO(4) <g> T ) /SU(2) + = SO(3) ® Z 2 , (11) 

with three Goldstone modes. T denotes the time rever- 
sal symmetry, and type B phase spontaneously breaks T . 
Therefore the GSM of phase B contains two disconnected 
sub-manifolds, with positive and negative Det[(Q)] re- 
spectively. 

The order of Q can be obtained through the follow- 
ing SO (4) invariant interacting Lagrangian for Dirac 
fermions on the honeycomb lattice: 

8 

c = J2(^M°-MQ t Q], (12) 

a=l 

The interaction — gtr[Q t Q] can be generated with 
SO (4) invariant interaction on the lattice, for instance 

S<i j> ' &j ~ — tr[<5*Q]/8 + • ■ • . A simple mean field 
calculation after the standard Hubbard-Stratonovich 
Transformation of Eq. [T^] shows that the type A phase 
has more favorable ground state energy compared with 
the type B phase on the honeycomb lattice. In the fol- 
lowing we will mainly focus on the analysis on the type 
A phase. 



III. PHASE DIAGRAM AROUND TYPE A 
PHASE 

Now we hope to understand the topological defects, 
and the phase transitions driven by topological defects 
in phase A. We will first give an argument about the 
phase diagram and phase transitions using the quantum 
numbers carried by the topological defects, and then a 
systematic description based on the Majorana liquid for- 
malism developed in ReL 3 will be presented. We will 
demonstrate that these two approaches match very well. 



A. Topological defects and phase transitions 

In phase A, since the GSM is [S 2 ® S 2 ]/Z 2 , both spin 
and charge sectors can have Skyrmion like defects charac- 



terized by homotopy group ^[S 12 ]. Again, let us assume 
that Q33 is the only component that acquires a nonzero 
expectation value, then (Q33) breaks the SO (4) symme- 
try down to residual symmetries generated by S z and T z 
in Eq.[5l According to Ref£~— , under our current assump- 
tion that N c || z (Qab is the QSH vector), a Skyrmion of 
the spin sector manifold S 2 p - ln carries charge 2e; and a 
Skyrmion current is identified as the charge current: 

J" = [ d 2 xe abc N:d„N b s d p N c s . (13) 

oir J 

For the same reason, a Skyrmion of the charge sector 
manifold will carry spin-1: S z = 1. For a general state 
with (Qab) 7^ 0, the spin-Skyrmion carries the quantum 
number of the U(l) c h ar ge residual symmetry, while the 
charge- Skyrmion carries U(l) sp i n quantum number i.e. 
spin and charge view each other as topological defects. 

The condensation of Skyrmions with nontrivial quan- 
tum numbers can lead to unconventional quantum phase 
transitions. For instance the proposal of deconfined 
quantum criticality is based on the observation that the 
Skyrmion of the Neel order carries lattice momentum^r— , 
hence the condensate of the Skyrmion is equivalent to 
the valence bond solid state. In our current case, since 
a charge-Skyrmion carries spin-1, if the charge- Skyrmion 
is condensed, then the SU(2) c h ar ge symmetry is fully re- 
stored, which implies that the condensate of the charge 
Skyrmion is a Mott insulator. Meanwhile, the residual 
spin symmetry is further spontaneously broken down to 
Z 2 <E> Z 2 . One of these Z 2 corresponds to changing the 
sign of electron, the other one corresponds to reversing 
the direction of N s . The GSM of the charge Skyrmion 
condensate is 

GSM = SU(2) spin /[Z 2 ® Z 2 \ - SO(3)/Z 2 . (14) 

We mod Z 2 from SO (3), because after the proliferation 
of charge Skyrmion, N c is completely disordered, and N s 
becomes a headless vector, due to the Z 2 in Eq. [TO] 

How do we determine the order of the charge-Skyrmion 
condensate unambiguously? As was pointed out in 
Refii 2 -, the phase of the 0(3) Skyrmion condensate can be 
identified as order parameters that share an 0(5) Wess- 
Zumino-Witten (WZW) term with the 0(3) order pa- 
rameter. Therefore, to unambiguously identify the order 
of spin-Skyrmion condensate, we need to seek for order 
parameters that have an 0(5) WZW term with vector 
¥ a = Qab- It turns out that the Neel order parameter 
N is the only candidate of the charge-Skyrmion conden- 
sate. For arbitrary 6, we obtain the following WZW term 
between <p a = Q a b and N ~ XP^^X- 

a=l " 

3i f 

- — / dud 3 xe abcde <t) a d x (t) b d v (j) c d T (j) d d u (j) e , 
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^ = ^ = xT a S b X , a =1,2,3, 

5 = 7V rf ~ XM^X, (15) 

Therefore the charge- Skyrmion condensate contains both 
headless vector N s and Neel order N, and N s _L N. Phys- 
ically the headless vector N s corresponds to spin nematic 
order S ab = 3N^N b - 6 ab (N s ) 2 , which is invariant under 
reversing the direction of N s . All these results will be 
confirmed later with the Majorana liquid formalism. 

Notice that manifold SO(3) is equivalent to the pro- 
jected manifold S 3 /Z 2 , which gives us a convenient way 
of parameterizing SO(3). Let us introduce SU(2) spinon 
z a with constraint \zi\ 2 + |z 2 | 2 = 1. This constraint im- 
plies that the SU(2) spinon z a parametrizes S 3 . Then by 
coupling z a to a Z 2 gauge field, the gauge invariant GSM 
of the condensate of z a automatically becomes SO(3) 13 . 
SO(3) manifold can also be viewed as the manifold of all 
the configurations of three perpendicular unit vectors qi , 
q 2 and 33. These three vectors can be parametrized as 

qj = z' a b z, q 2 + icfe = z t ia v a b z, (16) 

which automatically guarantees the perpendicularity of 
these vectors. In our situation, the three perpendicular 
vectors that characterize the GSM are N s , N and N s x N. 
Since N s is headless, the GSM is in fact SO(3)/Z 2 . And 
in the next subsection we will demonstrate that it is most 
convenient to describe this GSM by introducing a Z 2 ® Z 2 
or Z4 gauge field. 

Manifold SO(3) has homotopy group 7Ti[SO(3)] = Z 2 , 
therefore phase B has topologically stable half vortex. 
Using the CP(1) spinon description introduced in Eg. [TBI 
this half- vortex can also be viewed as the vison (a dynam- 
ical 7r-flux) of the Z 2 gauge field coupled to z a . Pictori- 
ally, a vison can be viewed as a configuration with (for 
instance) q\ being uniform in space, while q 2 and 9*3 have 
a vortex. Now since N s and N s x N are both headless 
vectors, this state also supports "half vison", where N is 
uniform, while N s has a half vortex in space. In fact, this 
half vison has a counterpart in phase A. Since phase A 
has GSM [S 2 <g)S 2 ]/Z 2 , there exists "double half vortex", 
and both N s and N c reverse direction after encircling this 
double half vortex. After phase A is destroyed by pro- 
liferating the charge-Skyrmion, this double half vortex 
becomes the half vison of SO(3)/Z 2 . 

In the spin-charge dual side of the theory, all the con- 
clusions can be obtained by straightforward generaliza- 
tion. Once the spin-Skyrmion is condensed, the sys- 
tem will also enter a phase with GSM SO(3)/Z 2 , and 
the SU(2) S pi n symmetry is fully restored, which implies 
that the condensate of the spin-Skyrmion is spin singlet. 
In Ref£, the authors proposed that after the prolifer- 
ation of Skyrmions of the QSH vector, the system en- 
ters a spin singlet s— wave superconductor. In our sit- 



uation, since there is a generic SU(2) c har g e symmetry, 
the s— wave superconductor is promoted to a phase with 
GSM SO(3)/Z 2 . If N c || z {Q ab is the QSH vector), the 
SO(3)/Z 2 manifold is characterized with headless vec- 
tor N c and s— wave superconductor. The order of head- 
less vector N c implies that the degeneracy between CDW 
and s— wave superconductor is spontaneously lifted. In 
general, the order after spin-Skyrmion proliferation can 
also be determined with the same WZW term analysis in 
Eq. [15] A full list of order parameters with WZW terms 
can be found in RefJ^. 

This Skyrmion condensation transition is described by 
the same CP(1) field theory as the deconfined quantum 
criticality 10 ™. How do we see the CP(1) transition di- 
rectly? The CP(1) model C = ^{(d^-iA^z] 2 describes a 
transition between a condensate of spinon z a and a pho- 
ton phase. The spinon condensate has GSM S 2 , while 
the photon phase has GSM S , as it is a condensate of 
the U(l) gauge flux. In our case the Skyrmion condensa- 
tion is a transition between [S 2 £3 S 2 ]/Z 2 and SO(3)/Z 2 , 
while SO(3) can be roughly viewed as S 2 <g> S 1 , therefore 
effectively the Skyrmion condensation is more or less also 
a transition between S 2 and S" 1 , so it is equivalent to the 
CP(1) transition. This hand- waving argument will be 
made precise in the next subsection by the Majorana liq- 
uid formalism. 

Since a spin-Skyrmion (charge-Skyrmion) carries 
charge-2e (spin-1), then the corresponding half-Skyrmion 
(vortex) will carry charge-e and spin- 1/2 respectively. If 
a charge-e excitation encircles around a spin- 1/2 excita- 
tion bound with a charge-vortex, the charge-e excitation 
will acquire a 7r phase shift; on the other hand, if a spin- 
1/2 excitation encircles a charge-e excitation bound with 
a spin- vortex, the spin-1/2 excitation will also gain a n 
Berry phase. This implies that in phase A charge-e and 
spin-1/2 excitations have mutual semion statistics. 

In phase A, a vortex is not a local excitation, and the 
gapless Goldstone mode of phase A makes the adiabatic 
braiding between two excitations impossible, therefore 
the semion statistics in phase A is not well defined. How- 
ever, later we will see that phase A is adjacent to a liq- 
uid phase where the spin-charge mutual semion statistics 
persists, and it becomes a well defined property. 



B. Phase diagram with Majorana liquid formalism 

From now on we hope to understand the phase dia- 
grams discussed above with a more solid formalism. In 
Ref4, we discussed a fractionalized phase of electrons by 
decomposing £ as following: 

C = z s z c x, 

Z s = <t> s + i<f>lS x + i<f> s 2 S v + i<i> s 3 S z , 

Z c = cj) c + i<j>lT x + i4> c 2 T» + i(b c 3 T z . (17) 
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FIG. 2: Phase diagram around type A phase with (Q33) 7^ 0. 



The electron £ decomposes into the bosonic fields Z s and 
Z c carrying its spin and charge respectively, and into 
the Majorana fermion X carrying the Fermi statistics. 
The resulting theory has a SO(4) s = SU(2) S>9 ® SU(2) c , g 
gauge invariance: Z s and X carry SU(2) Si9 charges, and 
Z c and x carry SU(2) Ci9 charges. 

After the operator decomposition, when both Z s and 
Z c are gapped out, one obtains the parent state, i.e. the 
Algebraic Majorana Liquid (AML) state with Lagrangian 



■Ad qa 



(18) 



The fractionalized Majorana fermion x nns the same 
mean field band structure as the physical Majorana 
fermion (. The gauge field A" and A a c li also couple 
to the spin and charge SU(2) rotors Z s and Z c as well. 
Since x 110 longer carries physical spin and charge quan- 
tum numbers, the fermion bilinears of x can oru y break 
the gauge symmetry, but not physical symmetry. If Z s or 
Z c condense, the formalism reduces to the two standard 
slave particle formalisms, with fermionic 2? or bosonic 
spinons^ respectively. 

Now let us assume x enters a type A phase i.e. the 
matrix field Q a b — xA a bX condenses. For instance let us 
take 



(Q33) = (x° z x) * 0. 



(19) 



Although the fractionalized Majorana fermion x nns the 
same mean field band structure as £, unlike the physi- 
cal QSH vector, nonzero {Qzb) breaks no discrete sym- 
metries (time-reversal, refection) when rotor fields Z s 
and Z c are gapped. This is because the gauge sym- 
metry released from gapping out the rotor fields can al- 
ways reverse the sign of (Q3&). This condensate of Q33 
breaks the SU(2) S ® SU(2) C gauge invariance down to 
U(l) <8> U(l) gauge symmetry generated by S z and 
T z . Sometimes it will be convenient to use the following 



4, 4Y 



'2 1 



(20) 



It was discussed in our previous work^ that, after inte- 
grating out x, we obtain the low energy theory when 
(Q33) 7^ 0, which is a mutual Chern-Simons theory: 

r — —p A z Pi a z 

~ 2tt up ^ p 



~s |2 



+ |(^-a*>-r + r.|«'| 

+ \(d^-iA z c ^z c a \ 2 + r c \z c a \ 2 



(21) 



The CP(1) fields z s a and z c a carry spin and charge respec- 
tively, and we have chosen the notation to make both 
spin and charge SU(2) physical global symmetries mani- 
fest. Eq. 1211 implies that there is a mutual semion statis- 
tics between the charge and spin CP(1) fields z c a and z s a , 
which verifies the observation in section 1111 Al 

This field theory is similar to the one on the triangular 
lattice^, with mutual semion statistics between charge 
and vison, except there the SU(2) symmetry of the vi- 
son is broken down to discrete symmetry by higher order 
terms^i, while here the SU(2) charge symmetry is exact. 
Based on this analogy, we can propose a similar global 
phase diagram with tuning parameters r s and r c (Fig. ^ 
as Refji^, with a different interpretation of the phases: 



Phase A 



Phase A is the phase with both z s a and z c a condensed, 
and the SU(2) s and SU(2) c gauge fields are both Hig- 
gsed and gapped out from the spectrum. This phase is 
characterized by the SU(2) vectors: 



N? 



N°: 



z s ^a a z s 



z c ^a a z c 



tv[ZlS a Z s S z ]- 
tr[ZlT a Z c T z ], 



(22) 



and the gauge invariant physical order parameter is 
Qab = £A ab C,=xZ t s Z t c A ab Z s Z c x 

~ tT[ZlT a Z c Tf} tr{ZlS b Z s S 9 } ( X A fgX ) 



(Q 33 )KN b c , 



(23) 



therefore the physical GSM is [S 2 ® S 2 ]/Z 2 . This phase 
is precisely the phase A obtained in the GL formalism 
in section [TTJ The Skyrmion of vector N s is equivalent 
to the flux of gauge field A z , and due to the mutual 
CS interaction, the gauge flux of A z s (Skyrmion of N s ) 
carries charge excitation , which confirms our analysis 
in section UlI Al 
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2. Phase A2 and A4 

Phase A2 has r s < 0, r c > 0, hence it is a phase with 
z s a condensed while z c a gapped. In section ITlI Al we con- 
cluded that the GSM of this phase is SO(3)/Z2, and it has 
both nematic order and Neel order by directly calculat- 
ing the WZW term. How do we understand the physical 
orders using the Majorana liquid formalism? Since z c a is 
gapped, in phase A2 there is no charge degrees of free- 
dom, therefore the "QSH" vector should correspond 
to a pure spin operator. In fact, since z c a is gapped, (Q a b) 
in Eq. [33] vanishes, hence the only gauge invariant oper- 
ator which acquires a nonzero expectation value is 

3Ql d Q de tr[ZlS a Z s S c ] tr[Z*S fc Z s S e ] 

~ (Q33> a JV?JV* (24) 

therefore we can define the physical order parameter as 

S ab ~ 3N°N b s - S ab (N s ) 2 ~ 3SyS* - 5 ab Si ■ S . (25) 

Hence S a b is the spin- 2 nematic order parameter which 
breaks the spin rotation symmetry down to U(l) ® Z2, 
but preserves the discrete symmetries. Notice that the 
physical order parameter is always a bilinear of N s . 

When z s a is condensed, the SU(2) s gauge field is Hig- 
gsed, then the SU(2) s gauge charge of x becomes equiv- 
alent to the physical spin quantum number of C after a 
SU(2) gauge transformation. If we take (Q33) ^ 0, the 
low energy field theory for fermions in phase A2 reads: 

C = Tpj^-iA^ip + miQ^-ija^, (26) 

where ip — \i + ^Xn- Based on the QSH physics, the flux 
of gauge field A\ carries spin: 

V X A z c ~ ^a*ij) tr[Z l s S a Z s S z }. (27) 

This equivalence between the flux and spin implies that 
the photon phase of gauge field A z c , which is the con- 
densate of the flux is a spin XY order. This effect was 
studied in Ref i 17 ' 18 with projected wave-function calcu- 
lation, and the photon phase of the U(l) gauge field is 
precisely the Neel order: 

N a ~ x/J, y S c x ti[ZlS a Z s S c ]. (28) 

Based on these analysis, we conclude that phase A2 is 
a phase with both nematic vector N s in Eq. [22] and AF 
Neel order N in Eq. [28] Eq. [22] and Eq. ES guarantee 
that N s ± N: 

N S -N~J2 ti[ZlS a Z s S z ] tr[ZlS a Z s S x ] = 0. (29) 

a 

As we mentioned before, since the nematic vector is head- 
less, the GSM should be SO(3)/Z 2 . This analysis again 
confirms our prediction in section IIII Al with the WZW 



term. 

Phase AA is the spin-charge dual phase of phase A2, the 
GSM is also SO(3)/i?2 with three branches of Goldstone 
modes. The spin-charge dual of the inplane Neel order 
is precisely a s— wave superconductor. The spin-charge 
dual of the nematic order S a b will break the SU(2) c h a rge, 
for instance 

S zz ~ 2/;/; - TfTf - T t y Tf, (30) 

with Ti given by Eq. [T] Therefore, if we turn on an extra 
density repulsion between next nearest neighbor sites in 
Eq. [T] it corresponds to turning on S zz , and breaks the 
SU(2) c har g e down to U(l) ®Z 2 - This U(l) corresponds to 
the ordinary electron charge conservation, and Zi corre- 
sponds to the discrete particle-hole symmetry. 

3. Phase A3 

Phase A3 is a liquid state with r s > 0, r c > 0, both z s a 
and z c a are gapped out. When z^ is gapped, A z c „ is in 
the photon phase. Since the photon phase of 2+ Id U(l) 
gauge field is also the condensate of gauge flux based 
on the standard QED-superfluid duality, the mutual CS 
coupling in Eq. [21] implies that the photon phase of A z c „ 
breaks the U(l) down to Z 2 gauge symmetry. For the 
same reason, U(l) is also broken down to Z^. The 
mutual CS theory in Eq. EHhas the same topological de- 
generacy as the standard Z 2 gauge field on the toru o 15 ' 19 , 
and the mutual statistics between charge and spin is an 
analogue of the mutual statistics between charge and vi- 
son of the well-known toric code models. 

In addition to the Z2 gauge field coming from the mu- 
tual CS coupling, there is one extra residual Z 2 gauge 
symmetry which corresponds to reversing the sign of 
gauge symmetry generators S z and T z simultaneously, 
while leaving Q33 invariant. This extra discrete Z2 gauge 
symmetry contains group elements 

G< 22 > =/ 4 x4, or S X T X . (31) 

This Z 2 gauge field couples to both spin and charge SU(2) 
rotors Z s and Z c , but it was not explicit in our continuum 
limit field theory. Therefore phase ^43 is characterized by 
Zi ® Zi gauge fields. Under Z 2 gauge symmetry G^ z2 \ 
the fractionalized particles transform as 

Z s ,j -> Z. d -((l + n J )I M + (l-ii j )iS*), 
Z cd -> Z cJ ^((l + Mj)i4x4 + (1 - Vj)iT x ), 

xj ^((l+^w + ti-iiimxi, 

Hj = ±1. (32) 

If the matter fields are ignored, these two Z 2 gauge fields 
are equivalent to a Z4 gauge field with group elements 
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G^ 4 ) = exp[i6S x T x ], 9 = 0, tt/2, tt, 3tt/2. 

These two Z 2 gauge field together again implies that 
the GSM of A2 (the condensate of z s a ) is SO(3)/Z 2 , as 
we already concluded. If we approach phase A2 from 
phase A3, we can interpret phase A2 as the condensate 
of SU(2) spin rotor Z s which couples to the two Z 2 gauge 
groups discussed above. With the condensate of Z s we 
can again define three perpendicular vectors: 

N s = ti[ZlSZ s S z ] ~ z s ^az s , 

Nj = tr[ZlSZ s S x ] ~ Re[(z s )Ha y az s ], 

N 2 = tr{ZlSZ s S y } ~ Im[(z s )Ha y az s ]. (33) 

N s and N 2 change sign under gauge transformation 
Z s — > Z s iS x (i.e. fij = —1), while Ni never changes 
sign, therefore N s and N 2 become headless nematic vec- 
tors by coupling to the Z 2 ® Z 2 gauge group. Hence the 
manifold formed with N s , N x and N 2 is SO(3)/Z 2 . This 
is completely consistent with our description of phase A2 
in section iHIAl and IIII B 2l 

In addition to the mutual statistics between Z s and 
Z c , there is one more topological defect in phase A3 with 

J[G {z2) = S X T X , (34) 
c 

C is a closed loop on the lattice. This defect car- 
ries a gauge flux S X T S . After encircling this defect, 
Z s -> Z s iS x , Z c -> Z c iT x . The vectors N s and N 2 
always acquire a minus sign after encircling this defect. 
Therefore this defect is a counterpart of the "double half 
vortex" and "half vison" discussed in section IIII Al 



4- Discussion 

The universality class of the phase transitions in Fig. [5] 
can also be analyzed with field theory Eq. in the same 
way as Ref Quoting the results in Ref the transi- 
tion between phases A3, A2, and the transition between 
phases A3, AA are 3d 0(4) transitions, because spinon z s a 
and z c a are 0(4) vectors, and the fully gapped discrete 
gauge fields coupled to the 0(4) vector do not change 
the 0(4) universality class^. The transition (A, A3), 
and the transition (A, AA) are CP(1) transitions, which 
become manifest with the CP(1) fields z s a and z c a , and 
U(l) gauge fields A z s ^, A%^. 

A recent quantum Monte carlo simulation of the Hub- 
bard model on the honeycomb lattice suggests that there 
is a fully gapped spin liquid phase 21 sandwiched between 
the ordinary Neel order and semi-metal phase, which has 
motivated spin liquid analysis on the honeycomb lattice 
using either slave boson or slave fermion technique o 22 ' 23 . 
In our formalism, phase A3 in phase diagram Fig. [5] is 
a candidate of this gapped spin liquid. However, based 
on our analysis, phase A3 is not directly adjacent to a 



pure Neel order, instead phase A3 is adjacent to phase 
A2 with both Neel order and nematic order. Starting 
with phase A2, we need to go through one more transi- 
tion which suppresses the nematic order, and enters the 
final Neel order in the large Hubbard U limit. 

The spin-2 nematic order is a natural candidate of 
the ground state of spin-1 systems, with bi-quadratic 
interactions 2 ^. For spin-1/2 models, nematic order can 
exist when there is a considerable ring exchange or multi- 
spin interaction, which can be generated in the weak 
Mott insulator phase of the simplest Hubbard model 
with high order perturbation of t/U. Our prediction 
of a phase with coexistence of nematic and Neel order 
can be checked numerically in future. We will present 
a general classification about nematic orders and their 
adjacent spin liquid phases in future 2 ^. 

A similar analysis can be applied to the Z 2 liquid phase 
obtained from the standard Schwinger boson formalism. 
Since the spinon z a always couples to a Z 2 gauge field, the 
condensate of z a is not an ordinary Neel order, because 
there always exists three perpendicular gauge invariant 
vectors like Eq. [TBJ In Ref^, the authors proposed that 
there is an intermediate chiral antiferromagnetic order 
between a fully gapped Z 2 liquid phase and a Neel order. 
This chiral AF state has GSM SO (3), which is different 
from the phase predicted in our paper with both nematic 
and Neel order. 



IV. PHASE DIAGRAM AROUND TYPE B 
PHASE 

Now let us move on the phase B with GSM SO(3)<g>Z 2 . 
All the phases discussed in this section breaks time- 
reversal symmetry T, therefore we only focus on one 
of the two disconnected sub- manifolds SO (3). Phases 
with GSM SO (3) have been studied extensively with non- 
collinear spin density wave^. After disordering the state, 
both SU(2) sp ; n and SU(2) c h ar ge are restored, while the vi- 
son of the SO (3) manifold is still locally conserved, and 
the system most naturally enters a Z 2 liquid phase. 

Again, we hope to understand the phase diagram with 
the Majorana liquid formalism. Let us assume \ m the 
parent state Eq. [15] enters the type B phase, for instance 
(Q11) = (Q22) = (Q33) 7^ 0. It would be convenient to 
introduce the following CP(1) fields 

z s = ( z f, z s 2 y = 

z c = {z\, z c 2 y = (0§ - i<f>i $1 - i^y. (35) 

The phase diagram around type B order is depicted in 
Fig. El 

Phase B in Fig. [3] with both z s a and z c a condensed is 
precisely the phase B in the Ginzburg-Landau descrip- 
tion in section HH with GSM S0(3) ® Z 2 . In phase B2, 
the SU(2) S gauge symmetry is Higgsed by the conden- 
sation of z s a , while the SU(2) C is broken down to Z 2 
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FIG. 3: Phase diagram around type B phase with (Qn) = 
(Q22} = <Q 33 > ? 0. 

gauge symmetry by (Q a b), this is the Z2 liquid phase we 
discussed above. The residual Z2 gauge symmetry is the 
subgroup of SU(2) c that changes the sign of X- Phase 
BA is the same Z2 liquid phase as phase B2. 

Now we turn to phase B3 in Fig. [5] In this phase 
both z 8 a and z% are gapped, while (Q a b) breaks the SO (4) 
gauge symmetry down to SU(2)+ = SU(2) s9 + SU(2) c>g , 
hence in this phase there is only one SU(2) gauge field 
A^G a : with G a = S a + T a . After integrating out the 
fermion \, the following Chern-Simons term is induced 
for the residual SU(2) + gauge field: 

C = — tr ( A A dA + -A A A A A ) 
Ait \ 3 / 

+ \(d,-Y,A«<T a )z S \ 2 +r s \z s \ 2 

a 
a 

A = A a ^a a dx^ : (36) 

therefore phase B3 is characterized by SU(2) CS theory 
at level 2, which is a nonabelian theory^. 

A different way of obtaining the same theory, is by 
turning on another order parameter XX m addition to 
Qab- The order parameter XX w ih drive the system into 
a quantum Hall state, and lead to the SU(2)i CS theory 
for both A" and A a c , which is similar to the CS effective 

theory of the chiral spin liquid stated. The order (Q ao ) 
requires A® = A® = therefore the final theory 
becomes the SU(2) CS theory at level 2 in Eq. GH Now 
by reducing the order (xx) to zero, the SU(2) CS theory 
is unchanged. More detailed properties of this phase will 
be further discussed in future^. 

In phase -B3, the residual gauge symmetry G Sjff is ei- 
ther G Ci9 or —G c .g. This Z2 structure does not show up 
in the Lie Algebra of the gauge group, but it implies that 
there is one extra Z 2 gauge field that couples to x an d 
either one of z s a or z c a . Hence when z s a or z c a condenses, 



the gauge field A a ^ is Higgsed, but the system still has 
a Z2 gauge symmetry, which characterizes the Z 2 liquid 
phase B2 and BA. 

The transition (B2, B3) and transition (B3, BA) is a 
Higgs transition, described by spinon z & a or z c a coupled 
with SU(2) CS theory in Eq. [36l The universality class 
of these transitions is not understood yet. 



V. SITUATION WITH SU(2) charge BROKEN TO 
U(l) 

charge 

When the SU(2) c h a rgo symmetry is broken down to 
U(l) c hargo ® Z2 symmetry which corresponds to charge 
conservation and particle-hole transformation, the degen- 
eracy between Q 3o and Q%b, Q2b is lifted. For instance, if 
an extra repulsive next nearest neighbor density interac- 
tion (linear with S zz in Eq. I30[) is turned on3^, the system 
favors to develop Q^b the system only has QSH order 
with GSM S 2 . Then according to Ref£, the Skyrmion of 
the QSH vector carries charge-2e, and Skyrmion conden- 
sate is a s— wave superconductor. 

If the system favors to have T-SC Qu, Qnb rather than 
Qsb, then depending on the microscopic parameters the 
T-SC can have orders with either Qn, || Q26 (type ^4) or 
Qu -L Q2b (type B). The type A phase has fully gapped 
fermion spectrum, with GSM~ [S 2 <g> S 1 ]/Z 2 . Here S 2 
corresponds to the spin direction of the triplet Cooper 
pair, while S 1 corresponds to the pairing phase angle. 
Again, there are spin and charge topological defects. For 
instance, the charge vortex (defect of the S 1 part of the 
GSM) carries spin-1/2 quantum number (quantum num- 
ber of the residual U(l) sp i n symmetry). The proliferation 
of the charge vortex leads to the phase A2 in phase di- 
agram Fig. [21 and the transition is a CP(1) theory with 
easy plane anisotropy on charge CP(1) field z c a intro- 
duced in Eq. 1201 which is consistent with the conclusion 
in Ref3. 

The proliferation of the spin-Skyrmion restores the 
SU(2) spin symmetry, but the U(l) c h ar ge symmetry is still 
broken. However, since the spin and charge sectors 
can change sign simultaneously without modifying the 
ground state, after the proliferation of the spin-Skyrmion 
there is still a Z2 gauge symmetry for the charge mani- 
fold. Therefore the GSM of this phase is S /Z2, which 
corresponds to a charge-4e superconductor, instead of a 
charge-2e superconductor. A similar scenario was dis- 
cussed for the polar state of the ultra-cold spin-f spinor 
condensate^, which also has the GSM [S 2 <g) S 1 ]/Z 2 . 

Unlike type A order, the type B phase with Qib _L 
Q2b does not have a fully gapped fermion spectrum. For 
instance, with (Qu) = (Q22) 7^ 0, only spin-up is paired 
and gapped out, while spin-down is not gapped. Since 
the fermion spectrum is gapless, the quantum number of 
defects is no longer topologically stable. 
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VI. SUMMARY AND DISCUSSION 

In this work we have classified the QSH and T-SC 
states on the honeycomb lattice using the SO (4) sym- 
metry for a large class of extended Hubbard models. By 
analyzing the quantum numbers of topological defects, 
we obtained two different phase diagrams, which were 
also confirmed by the Majorana liquid formalism. The 
results of this paper can be straightforwardly generalized 
to other bipartite lattice. Our formalism also predicts 
a phase with both spin nematic and Neel order, sand- 
wiched between a fully gapped Z2 ® Zi liquid phase and 
the ordinary Neel order, which can be checked in future 
using the similar method as ReL^I. 

The results we obtained in this work explicitly demon- 
strates the spin-charge duality of the Hubbard model. 
For instance, in phase diagram Fig. [3J spin and charge 
view each other as topological defects. A similar spin- 
charge duality was applied to the cuprates high temper- 
ature superconducto r 31 ' 32 , and a global phase diagram 
with both spin and charge excitations was studied re- 
cently in Ref.— . We also note that other author o 18 ' 26 



have also studied the duality between spin and charge 
with the presence of QSH order parameters, for instance 
an easy plane version of spin-charge duality was identified 
as the self-duality of the easy plane noncompact CP(1) 
theory, in a model with inplane spin anisotropy. This 
duality led to a direct transition between inplane Neel 
order and d— wave superconductor. In our current work 
we showed that the generic symmetry of the Hubbard 
model and the condensate of matrix order parameter Q 
in Eq. [6] give us a complete and explicit duality between 
spin and charge in interacting electrons. 

In both Fig. [2] and Fig. [3] there is a multi-critical point 
with r s = r c = 0. The multi-critical point in Fig. [2] was 
analyzed in Ref.~^, and for large enough spinon number 
this multi-critical point is stable. Also, it has been pro- 
posed that a similar multi-critical point is responsible for 
the spin liquid behavior in material k — (ET) 2 Cu 2 (CN) 3 
on the triangular lattice 3 -^. The multi-critical point in 
Fig. [3] is more complicated, we will leave this multi- 
critical point to future study. 

The author appreciates the very helpful discussions 
with Leon Balents and Andreas Ludwig. 
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